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Abstract. We review recent progress on lattice calculations of nucleon spin structure, including the 
parton distribution functions, form factors, generalization parton distributions, and recent develop- 
ments in lattice techniques. 

Keywords: Lattice QCD, nucleon spin structures, form factors, GPDs 
PACS: 12.38.-t 12.38.Gc 13.40.Gp 14.20.Dh 13.60.Fz 

INTRODUCTION 

Quantum chromodynamics (QCD) has been successful in describing many properties 
of the strong interaction. In the weak-coupling regime, we can rely on perturbation 
theory to work out the path integral which describes physical observables of interest. 
However, for long distances perturbative QCD no longer converges. Instead, we use a 
discretization of space and time in a finite volume to calculate these quantities from first 
principles numerically; such research forms the regime of lattice QCD[[I1]. 

To keep the systematic error due to discretization under control, one follows Symanzik 
improvement order by order in terms of the ultraviolet cutoff (a) for both the action and 
operators. However, the breaking of continuous (Euclidean) SO(4) symmetry allows 
many new degrees of freedom, leading to various lattice actions that return to the same 
continuum action once the symmetry is restored. Thus, there exist many gauge and 
fermion actions for us to choose from. Today, most gauge actions used are 0(a 2 )- 
improved and leave small discretization effects (0(a 3 AQ CD )) due to gauge choices. 
On the other hand, most fermion actions are only 0(a)-improved and have systematic 
errors of 0(a 2 Ag CD ) that become dominant. For this reason, lattice calculations are 
generally distinguished according to the fermion action used. Differences among the 
actions are benign once all systematics are included, and the choice of fermion action 
is constrained by limits of computational and human power and by the main physics 
focus. The commonly used actions are: domain-wall fermions (DWF), overlap fermions, 
Wilson/clover fermions, twisted-Wilson fermions and staggered fermions. 

Since the real world is effectively continuous and infinitely large, we will have to take 
limits of a — > and V — > °° to eliminate the artifacts introduced in a discretized finite 
box. With the most state-of-the-art supercomputer, we are close but yet to simulate at 
the physical pion mass. Using calculations at multiple heavier pion masses, which are 
affordable for available computational resources, we can apply chiral perturbation theory 
to extrapolate quantities of interest to the physical limit. A recent work by the BMW 
collaboration [2] calculating multiple lattice spacings, volumes and pion masses as light 



as 180 MeV provided an excellent demonstration of how ground- state hadron masses 
with fully understood and controlled systematics are consistent with experiment. Such 
calculations with multiple pion masses also help to determine the low-energy constants 
of chiral effective theory. 

A typical nucleon interpolating field used in the lattice calculation is Xn = 
Y.x,a,b,c eip X£ [ U a^-Y5^b] u c , and the nucleon two- and three-point Green functions are 
obtained from 

r (2) (w) = (XN(t)xU^c)) (i) 

where & is the operator of interest. For the vector (axial) current, the operator is 
6 = ^7^(75) V- F° r the structure functions, the operators are 

(*% 

where D = l(D — D) is the difference between forward and backward co variant 
derivatives. We calculate only the "connected" diagrams, which means the inserted 
quark current is contracted with the valence quarks in the baryon interpolating fields, 
as in the majority of lattice three-point calculations. However, due to isospin symmetry, 
isovector quantities have a cancellation that removes the unknown disconnected piece. 
(Disconnected diagrams are notoriously difficult to calculate directly on the lattice. They 
require that expensive fermion-matrix inversion be applied to numerous source vectors. 
Much effort has been devoted to solving this difficulty in the near future with new 
techniques.) For more details, please refer to a selection of plenary talks: Refs. |30 and 
references within. 



LATTICE QCD REVEALS THE STRUCTURE OF THE NUCLEON 

The nucleon axial charge is well measured in neutron /3-decay experiments, so it is a 
natural candidate for demonstrating how well lattice QCD can be extrapolated to the 
physical pion point. The isovector axial charges gA are defined as the zero-momentum- 
transfer limits of the isovector axial form factors. Results from various collaborations 
are shown on the left-hand side of Figure [Q We show a small-scale-expansion fit 
(gray band) to the 2+1 -flavor (DWF sea quarks) RBC data[5] (blue filled circles). 
The results are consistent with the LHPC dataflg] (DWF valence with staggered sea); 
RBC's 2f and Of results are consistent with QCDSF's 2f Clover and Of overlap fermion 
numbers[|71] respectively. The lowest-pion-mass values among RBC's 2+lf and 2f results 
may suffer from sizable finite-volume systematic errors; larger- volume calculations 
should be carried out to confirm these suspicions. 

Results for the zeroth moment of transversity (l)§ q are given on the right-hand side 
of Figure [T] We observe rather weak dependence (roughly linear) on the quark mass, 
which remains consistent for calculations with the same number of sea-quark flavors. 
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FIGURE 1. (left) Renormalized axial charge versus pseudoscalar mass squared from various lattice 
groups, (right) Zeroth moment of transversity from different lattice groups. The bands in both plots are 
chiral extrapolations fit to the RBC 2+ If data. 



The chiral extrapolation form[4] is applied to the RBC 2+ If data[5], yielding 0.56(4) 
at the physical pion mass. However, this is close to what has been found by LHPC 
with mixed action[6]: about 0.7. These extrapolated values are significantly smaller than 
those found at the simulated pion masses, which are near 1. We urgently need data at 
the lightest pion mass to confirm the rapidly decreasing behavior predicted by the chiral 
effective theory. 

Figure |2] shows the latest calculations of the first moments of the quark momentum 
fraction (left) and helicity (right) distributions. Here again the Of, 2f and 2+lf results 
are compared between different collaborations with different choices of fermion actions 
and are seen to be consistent among calculations with the same number of sea-quark 
flavors. The chiral extrapolation Q4f| is performed using RBC's 2+lf data, which gives 
0.133(13) and 0.203(23) for the first moments of the quark momentum fraction and 
helicity distributions respectively, consistent with experiment. We see strong curvature 
due to the chiral form; more light-pion points should be taken to reduce extrapolation 
uncertainties. These extrapolation numbers are also consistent with the LHPC's mixed- 
action calculation^. 

Higher moments of the isovector distributions are calculated by LHPC (SCRI, 
SESAM) with Of and 2f Wilson and clover fermions[@] and by QCDSF with Of clover 
fermions at multiple lattice spacings[7]. Consistent results are seen among different 
groups: The second and third moments are about 25% and 10% of the first moments 
respectively. However, for moments with n > 4, divergences occur involving lower- 
dimension operators at finite a, which limits the number of moments accessible to lattice 
QCD. 

The twist-3 first moment of the polarized structure function d n is another interesting 
feature to consider. It is related to the polarized structure functions g\ and gj and the 
Wandzura-Wilczek relation d^. The lowest moment d\ from RBC's 2+lf data extrapo- 
lated to the physical pion mass is consistent with zero d\ aie = —0.002(2). Combined 
with the small value of J2 found by QCDSF's 2f calculation[7], we conclude that the 
Wandzura-Wilczek relation between the moments of gi and g2, which asserts vanishing 
d n , is at least approximately true. 

Studying the momentum-transfer (Q 2 ) dependence of the elastic electromagnetic form 
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FIGURE 2. Global comparison of the first moments of the quark momentum fraction (left) and helicity 
(right) distributions in terms of M\ and their chiral extrapolations. The bands in both plots are chiral 
extrapolations fit to the RBC 2+ If data. 



factors is important in understanding the structure of hadrons at different scales. There 
have been many experimental studies of these form factors on the nucleon. A recent such 
experiment, the Jefferson Lab double-polarization experiment (with both a polarized tar- 
get and longitudinally polarized beam) revealed a non-trivial momentum dependence for 
the ratio G P E / G». This contradicts results from the Rosenbluth separation method, which 
suggested \i p G E jG p M ~ 1. The contradiction has been attributed to systematic errors due 
to two-photon exchange that contaminate the Rosenbluth separation method more than 
the double-polarization. (For details and further references, see the recent review arti- 
cles: Refs. 09D.) Lattice calculations can make valuable contributions to the study of nu- 
cleon form factors, since they allow access to both the pion-mass and momentum depen- 
dence of such form factors. Recently, the limitations of the largest-available Q 2 (in terms 
of the quality of the signal-to-noise ratios) has been overcome ||10|l . An exploratory study 
using clover fermions extends the range of momentum transfer to 6 GeV 2 , as shown in 
Figure |3l The Q 2 dependence of the neutron has exceeded the range of the current exist- 
ing data. Such calculations will provide interesting comparisons for data collected after 
the future 12-GeV upgrade at Jefferson Lab. 

The generalized parton distributions (GPDs) have been calculated by LHPC (2+ If 
mixed action, M n ~ 350-760 MeV)|@] and QCDSF (2f clover action, M K ~ 340- 
950 MeV)[7]. One of the main topics of physics interest derived from GPDs is to study 
the origins of the nucleon's spin. The quark spin (AE 9 ) and total quark contribution (J q ) 
to the angular momentum can be connected to GPDs via [jll|l 

\^ = lAJ (0); J« = ^[Af (0) +4,(0)], (3) 

where A q nQ and {A,B} q nQ are the polarized and unpolarized generalized form factors. The 
left-hand side of Figure 0] shows the quark spin and orbital angular momentum (L q ) 
calculated by LHPC and QCDSF; they are consistent. One found the total d quark 
angular momentum and the angular momentum from the sum of u and d quarks to 
be consistent with zero. This is because L d and L u are of the same magnitude but 
have opposite signs. A similar relation holds for AZ d and L d . Note that in both of the 
calculations only connected diagrams are included. 
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FIGURE 3. Nucleon form factors with pion masses of 480 (triangles), 720 (circles) and 1080 (dia- 
monds) MeV. The dashed lines are plotted using experimental form-factor fit parameters^. 



For the transverse structure of the proton, QCDSF calculates the lowest two moments 
of the transverse spin densities of quarks in the nucleon[7]. An exploratory attempt 
has been made to calculate transverse-momentum distributions (TMDs) using lattice 
link products to approximate the Wilson line in light-cone frame in spatial coordinates 
and Fourier transforming into momentum-space to map the density distributions[@]; this 
calculation used 2+lf mixed action (DWF on staggered, M n ~ 500 MeV). 

Another new development during the last year is progress on disconnected diagrams. 
An indirect way of calculating the strangeness form factors via charge symmetry in 2+lf 
lattices is shown on the right-hand side of Figure |4] The direct approach has also made 
great progress. For example, the strange-quark distribution (x) s has been calculated by 
the ^QCD collaboration on 2+lf lattices[12]. The same group has also calculated the 
gluon momentum fraction (x) g with improved signal. 

This is an exciting era for the use of lattice QCD in nuclear physics: there have been 
huge leaps due to increasing computational resources worldwide and improved algo- 
rithms, allowing continual improvement in lighter pion masses, larger volumes and finer 
lattice spacings. Various groups have demonstrated universality with the consistent re- 
sults coming from independent calculations using different lattice actions. By reproduc- 
ing well measured experimental values, we solidify confidence in lattice predictions of 
quantities that have not or cannot be measured by experiment. There are many different 
aspects of hadron structure that one can do with lattice QCD; only a few examples have 
been presented here. 

In the near future, pion masses around 200 MeV or lighter with multiple volumes 
and lattice spacings will become commonplace. There will be less need to depend on 
chiral perturbation theory for extrapolations, and once the physical pion mass becomes 
accessible, we can check its correctness. Full-contraction calculations (including discon- 
nected diagrams) in all matrix elements, form factors and GPDs will lead to precision 
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FIGURE 4. (left) The quark-component contributions to quark spin and orbital angular momentum 
in the spin of the nucleon from LHPC (squares, diamonds and upward triangles) and QCDSF (tilted 
squares and rightward triangles). The filled/open symbols represent dynamical/quenched data. The bands 
are chiral extrapolations based on LPHC's dynamical data, (right) Experimental results (ellipses) and 
various theoretical calculations (points) of G S M -G S E 



calculations for individual quark components or individual proton and neutron quanti- 
ties. Further improvements in methodology and expanding computational resources will 
allow direct calculations of the gluon content of the nucleon. The future is unlimited 
with lattice QCD. 
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